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Abstract. In this paper we use a probabilistic approach to derive the expres- 
sions for the characteristic functions of basic statistics defined on permutation 
tableaux. Since our expressions are exact, we can identify the distributions of 
basic statistics (like the number of unrestricted rows, the number of rows, and 
the number of Is in the first row) exactly. In all three cases the distributions 
are known to be asymptotically normal after a suitable normalization. We also 
establish the asymptotic normality of the number of superfluous Is. The latter 
results relies on a bijection between permutation tableaux and permutations 
and on a rather general sufficient condition for the central limit theorem for 
the sums of random variables in terms of dependency graph of the summands. 



1. Introduction 

Permutation tableaux are relatively new objects that are in bijection with per- 
mutations |Bul ICoNl [SW . They were introduced in the context of enumeration 
of the totally positive Grassmannian cells [0 IWj . More recently, permutation 
tableaux generated additional research activity when they have been connected in 
[Col ICoWl ICoWlj to a particle model in statistical physics called the Partially 
ASymmetric Exclusion Process (PASEP); see [BrEl IBrCPREl IDDMl [DEHP, 
IDSI [Sj IUSW] for more information on PASEP. 

A permutation tableau [SW] is a Ferrers diagram of a partition of a positive 
integer into non-negative parts whose boxes are filled with Os and Is according to 
the following rules: 

(1) Each column of the diagram contains at least one 1. 

(2) There is no which has a 1 above it in the same column and a 1 to its 
left in the same row. 

An example is given in Figure [T] 

The size parameter of a permutation tableau is its length defined as the number 
of rows plus the number of columns of the tableau. For example, the tableau in 
Figure [T] has 6 rows and 7 columns so its length is 13. 
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Figure 1 . Example of a permutation tableau 

Different statistics on permutation tableaux were denned in [CoWll ISW) . We 
recall that a in a permutation tableau is restricted if there is a 1 above it in the 
same column. A row is unrestricted if it does not contain a restricted 0. A 1 is 
superfluous if it has a 1 above itself in the same column. We will be interested in the 
number of unrestricted rows, the number of superfluous Is, as well as the number 
of Is in the top row, and the number of rows. For example the tableau in Figure 
[T]has three superfluous Is, six rows, three Is in the top row, and five unrestricted 
rows. 

Most of the past research on statistics of permutation tableaux was based on 
bijections between permutation tableaux of length n and permutations of [n] := 
{1, . . . , n} and using known properties of permutations. This was not always easy 
as sometimes it is not that easy to see into what parameter a given statistic is 
mapped by a bijection. In [CoHj a direct approach based on a probabilistic con- 
sideration was proposed. It enabled the authors to compute the expected values of 
these statistics in a simple and unified way. In this work we go one step further and 
we will concentrate on the limiting distributions of (properly normalized) statistics 
mentioned above. In most cases we will accomplish it by computing the probability 
generating functions of the quantities in question. As a matter of fact, since the 
approach proposed in [CoH allows for exact (and not only asymptotic) computa- 
tion, in most cases we will be able to identify the distribution of a given statistics 
exactly and not only asymptotically. As a consequence, we can see, in particular, 
that: 

• the number of unrestricted rows in a random permutation tableaux of 
length n has the same distribution as the number of cycles, or the number 
of records, in a random permutation of [n] (a record in a permutation 
a = ((Tj) is any <Tj such that o~i > <Tj for j < i); 

• the number of rows has the same distribution as the number of descents 
in a random permutation counted by Eulerian numbers (a descent in (pi) 
is any pair (<jj, <7i+i) such that er,; > <7j + i); 

• the number of Is in the first row is, in distribution, one less than the 
number of unrestricted rows. 

The first two results follow from the analysis of known bijections between permuta- 
tions and permutation tableaux (see [CoNllSW] ). while the third is a consequence 
of an involution on permutation tableaux presented in [CoWlj . However, neither 
the involution from [CoWlj nor any of the bijections from [CoNl ISWj are very 
straightforward. In addition, the first two results given above required two different 
bijections in [SWj . So, while these bijections carry more information, they turn 
out to be rather cumbersome to work with. Looked at from this perspective, our 
approach provides a streamlined and unified derivation of the above results. 
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As for the number of superfluous Is, the situation is a bit more complicated. 
Although we do derive its probability generating function, we do not deduce directly 
from this the asymptotic normality of its distribution. Instead, in this case we 
will rely on the fact that one of the two bijections between permutation tableaux 
and permutations described in [SW] sends the number of superfluous Is to the 
number of generalized patterns 31—2 and we will prove the central limit theorem 
for the number of such patterns in random permutation of [n] (an occurrence of a 
generalized pattern 31—2 in (at) is any pair 1 < i < j such that > Uj > 

<7i). Our proof is based on a rather general sufficient condition for the central limit 
theorem developed in [J] and [JLR1 Section 6.1]. Bona |Bo| has used the same 
method to show similar results on asymptotic normality of general permutation 
patterns, but as far as we know the central limit theorem for the number of a 
generalized pattern such as 31—2 is new. Furthermore, the number of permutation 
tableaux of length n with k superfluous Is is equal to the number of permutations 
of [n] with k crossings (a crossing |Coj in a permutation (di) is a pair such 
that i < j < Cj < <jj or i > j > Oj > tjj). It follows, therefore, that the number of 
crossings in a random permutation of [n] is asymptotically normal. 

2. Basic facts 

Let T n be the set of all permutation tableaux of length n. We denote the uni- 
form probability measure on T n by P„ and E n will denote the expectation with re- 
spect to P„. We denote by R ni C n , U n , F n , and S n the random variables represent- 
ing the numbers of rows, columns, unrestricted rows, Is in the top row, and superflu- 
ous Is, respectively, in a random tableau of length n. (We may allow n = 0; there is 
a single empty permutation tableau of length 0, and Uq = Fq = Rq — Cq = Sq = 0. 
It is sometimes convenient to start inductions with n = 0, but for simplicity we 
often treat n > 1 only.) 

If is a er-algebra then E( • \T) denotes the conditional expectation given T . 
If X is non-negative integer valued, for example one of the permutation tableaux 
statistics just defined, we let gx{z) = Ez x be its probability generating function 
(in general defined at least for \z\ < 1; for the variables consider here, z can be any 
complex number). We will often omit the subscript X. 

The arguments in |CoH) were based on a construction of tableaux of size k 
from tableaux of size k — 1 by extending the length of the latter from its south-west 
(SW) corner either to the south (creating a new empty row) or to the west (creating 
and then filling a new column). Note that each permutation tableau in Tk is an 
extension of a unique permutation tableau in Fk-\- We refer the reader to |CoH[ 
Section 2] for details. We let Mk indicate the direction of the move at the kth step 
(i.e. when the length is increased from k— 1 to k). We refer to Mk as the fcth move 
and will write Mk = S or Mk — W to indicate its direction. 

The following simple observations were crucial for the arguments in |CoHj and 
are crucial here as well (see |CoH] for some details that are omitted here). 

• When extending a tableau by a W move, the new column has to be filled 
with 0s and Is; each restricted row must be filled with 0, but the unre- 
stricted rows can be filled arbitrarily except that there must be at least 
one 1. Hence, a tableau with U unrestricted rows can be extended by a 
W move in 2 U — 1 ways; since there always is a unique S move, the total 
number of extensions is 2 U . 
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• We use T n to denote a generic element of T n , and use T n -i to denote 
the corresponding element of T n -i (i.e., such that T n is an extension of 
T n —i); we let U n = U n (T n ) be the number of unrestricted rows in T n 
and U n -i = [/„_i(T„-i) the number of unrestricted rows in T n _i, and 
similarly for the other variables that we study. In this way, P„ induces a 
probability measure (also denoted P n ) on T n -\, namely each element of 
T n —i is assigned a measure that is proportional to the number of tableaux 
from T n it generates when its length is increased from n— 1 to n. Note that 
this differs from the uniform distribution P„_i on T n —i. Since \T n \ = nl 
and a tableau T from 7^_i generates 2 Un ~ 1 ^ tableaux of length n the 
relationship between these two measures is easy to find. One way to state 
this relationship is that if X is any random variable on T n -\ then 

(2.1) E n X = -E n ^(2 u "-iX) , 

n 

where integration on the left is with respect to the measure induced on 
7n-i by P„ and integration on the right is with respect to P n -i- Equiv- 
alently, if we let J- n -i be the er-algebra on T n generated by the mapping 
T n h-> T n _i, then dP„/dP„_i = 2 u "^/n on T n - X . 

• The sequence of the distributions of the number of unrestricted rows is 
given by Uq = 0, U\ = 1 and for n > 1 the conditional distribution of U n 
given T n -i (i.e., given T n -i) is (under P„) 

(2.2) £{U n | T n - X ) = 1 + Bin(t^_i), 

where Bin(m) denotes a binomial random variable with parameters m and 
1/2. 

To illustrate how these facts are put together to work consider the number of 
unrestricted rows, U n . For its probability generating function we have: 

g Un (z) = E n z u " - E„E„ (z u ~ \ ^ n _ x ) = E„E„ ^i+Bin(t/„_0 | Tn 



(2.3) =zE n — — =-E„ 

\ 2 / n 

= -E„_ 1 (z + l) C/ - 1 , 
n 

where we have used (in that order) conditioning, (|2.2p . the obvious fact that 
SBin(m)(*0 = Ez Bin ( m ' = ( £ ^ i )'™, and (|2~T|) . It follows by induction, and guAz) = z, 
that 

, , T(z + n) I z + n\ '-J-r 1 z + j -A- / 1 z 

^ ^^=w = ( « J = njTi=n(i-^+^ 

W V 7 j=0 J fe=l 

The factor on the right-hand side is the probability generating function of a random 
variable that is 1 with probability l/k and with probability 1 — 1/k. Since the 
product of probability generating functions corresponds to summing independent 
random variables, we obtain the following statement 

Theorem 2.1. The number of unrestricted rows U n is distributed like 

n 

fe=i 
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where 3\, J2, ■ • ■ , are independent indicators with P(Jk = 1) = l/k. In particular, if 

h n = Y^k=i ^ 1 an d ^»* 2 ' = Sfc=i ^ _ 2 * s ^ e harmonic and generalized harmonic 
number, respectively, then 

n n n 1 / 1 \ 

E n [/ n = ^P(</ fe ) = /i„, var„(C/„) = ^var(J fc ) = ^- (i-t) =h n -h&, 

k=l k=l k=l * ' 

and 

^i,JV(0,l), 
v Inn 

where N(0, 1) denotes a standard normal random variable. 

Remark 2.2. It is seen from the above statement that the distribution of U n 
coincides with the distribution of the number of cycles in a random permutation 
of [n] (see for example [FJ Chapter X, Section 6(b)] or [R] Chapter 4, Section 3]), 
or, equivalently, the number of RL minima in a random permutation (which is 
known to be equidistributcd with the number of cycles). Indeed, a bijection between 
permutation tableaux and permutations described in [CoN maps unrestricted rows 
onto RL minima in the corresponding permutation. 

Remark 2.3. More generally, we may define on T n permutation tableaux T% £ 
Tk for every k < n, such that Tk+\ is an extension of and let Tk be the er-algcbra 
generated by Tk- It can be seen by induction, similar to (|2.3[) . that a permutation 
tableau Tk of length k can be extended to (n — k) \ (n — k + 1) C/ »=( T '=) permutation 
tableaux in T ni and hence dP n /dPk — (n — k + l) c/fc /(^) on Tk : 1 < k < n. 
Furthermore, (|2.4p can be generalized to 

En(^" I Fk) = r(Z + V, fc)fc! (^ + n kr ■ 
1 (z)n! 

We can further study the sequence (Uk)k = i as a stochastic process under P„; similar 
calculations show that this is an inhomogeneous Markov process with transitions 
given by C(U k+1 | «F fc ) = 1 + Bin(L4, ra ^ x ), which generalizes 



3. The number of Is in the first row 

Recall that F n denotes the number of Is in the top row of a random permu- 
tation tableau of length n. The following fact is a consequence of an involution 
on permutation tableaux given in [CoWlj . Since this involution is not too easy 
to describe we provide a more direct justification of (|3.1[) based on the arguments 
given in the preceding section. 

Theorem 3.1. For every n > 1, we have 

(3.1) F n = U n - 1. 

In particular, 

E n F n = h n -l, var„(F„) = h n ~ h£\ 

and 

F„ — In n a 
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To establish (|3.1|l it will be convenient to prove some auxiliary lemmas. Define 
Gk to be the position (counting only unrestricted rows) of the topmost 1 on the 
fcth move, provided that this move is W. Then 

n 
k=X 

Note that Gk is undefined if the fcth move is S. While it is inconsequential in 
this section as we will be interested in the event {Gk = 1}, for the purpose of the 
subsequent sections it is convenient to set Gk = Uk-i + 1 if the fcth move is S. We 
have 

Lemma 3.2. For all fc > 2, 

Pk(G k = 1 | T k -i) = \. 

Proof. A given tableau T e Tk-x has 2 Uk - 1 ^ extensions. The number 
of extensions by a W move with a 1 in the topmost row (which always is unre- 
stricted) is 2 Uk - 1 ( T )~~ 1 since the other Uk-i(T) — 1 unrestricted rows can be filled 
arbitrarily. Hence the probability that Gk = 1 in a random extension of T is 
2tf*-i(T}-i/ 2 tf»-iCr) =1/2. □ 

Lemma 3.3. For any complex numbers z,w and every k>2, 

c . /,• / i r , z + w {w + l^ tl ' 
E k {z Gk w Uk T k -i) = w — — 
w + 1 

PROOF. Let T e Th-\- If T is extended with Gk = 1, i.e., by a W move with 
1 added in the topmost row, then we may add or 1 arbitrarily to all other unre- 
stricted rows, and the rows with 1 added remain unrestricted. Hence, conditionally 
given Gk = 1, the extension has 1 + Bin(C/fc_i(T) — 1) unrestricted rows. Recall 
from (|2.2[) that without further conditioning, the number of unrestricted rows is 
1 + Bin({/fc_i(T)). Consequently, using also Lemma EOI 

E k (z Ia *w Uk | F k -i) = zE k {w Uk I Gk=1 | Tk-x) + E k {w u H Gk ±x \ Tk-x) 
= {z- l)E k {w u "I Gh=1 | Fk-i) + £ k {w Uk | T k -x) 
= (z- l)E k (w 1+Bia ^-^I Gk=1 | F k -x) + E k (w 1+Bi ^ u ^ | Fk-x) 



( z _l) w (__) P k (G k = 1 | Tk-x 



ii- 



io /'w + l\ Uk - 1 fw + l^ Uk - 1 



w + 1 V 2 

u k - 



w + 1 

Lemma 3.4. The joint probability generating function of F n and U n is given by 

r- i F tt \ T(z + w + n — 1) T-r z + w + k — 2 

E n {z "w u ") = w , , = w M . 

I (z + w)n\ ±A - k 

v ; k=2 
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Proof. By Lemma 1331 and (|27T|) . 



E, 



z + w 1 



In 



and the formula follows by induction. □ 
We now can complete the proof of Theorem 13.11 
Proof of Theorem 13. 11 Taking w = 1 in Lemma \3A\ we obtain 

p f„ _ T(z + n) 
nZ r(z + l)n!' 

which equals E n z Un ~ x by (23). □ 



Note that we recover (12.4[) by taking z = 1 in Lemma 13^41 We can also describe 
the joint distribution of F n and Z7 n . 

Theorem 3.5. For every n > 1, i/ie jom< distribution of F n and U n is given 

by 

n 

(F„,C/„) = X^' Jfc )' 
fc=i 

where the random vectors (Jk,Ik) are independent and J\ — 0, ii = 1, and, /or 
* > 2, P(/ fc = 1, J fe = 0) = P(J fe = 0, J fe = 1) = P(J fe = 0, J fe = 0) = 1 - 2/k. 

PROOF. We have E(z Jk w Ik ) — (k - 2 + z + w)/k if k > 2 and tt> if fc = 1, 
and thus the joint probability generating function of the right-hand side equals the 
product in Lemma 13.41 □ 

Corollary 3.6. The covariance of F n andU„ is 

cov n (F„, U n ) = E n (F n U n ) - E n F n E n U n = -(h^ - 1). 

PROOF. We have cov(J fe ,/ fc ) = — EJ fe E/ fe = — fc~ 2 for k > 2, and cov(Ji,ii) = 
0, and thus by Theorem 13.51 cov(F„ , U n ) = J2k=2 ^~ 2 - ^ 

It follows further easily from the central limit theorem for vector- valued random 
variables, e.g., with Lyapunov's condition, that the normal limits in Theorems 12.11 
and 13. II hold jointly, with the joint limit being a pair of two independent standard 
normal variables. We omit the details. 

4. The number of rows 

In this section we consider the number of rows R n in a random permutation 
tableau of length n. Since the rows of a permutation tableau correspond to S steps 
in the process of its construction we can write 



k=l 
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Lemma 4.1. For 1 < k < n and any complex numbers z, w we have 

Proof. By conditioning on T k -i we get 

(4.1) E k {z R *w u ») = E k (z^Ekiz^'^w^ | T k ^)) . 

Note that M k = S if and only if Uk = 1 + U k -i and, conditionally on T k -\, this 
happens with probability 1/2 Hence, using also Lemma [3731 with z = 1 (or the 
argument in (|2.3p ). 

Efc^ 1 "*-*^* | F k -i) = E fe ((z^= s - l)w u " | + E k (w u " | T k ^) 



= E k ((z-l)I Mk =sw Uk \Fk-i) 

x w 1+u *-i fw + l\ Uk - 
= (z - 1) 



2U k -i 1 

=^-Dy +w {—) • 

Putting this into (|4.ip and applying (|2.1j) proves the lemma. □ 
We can now compute the probability generating function of R n . 
Theorem 4.2. We have 

where (?) are i/ie Eulerian numbers counting the number of permutations a — 
(<ri, . . . , cr n ) o/ [n] uwi/i /c descents (recall that a descent is a pair (aj, ffj+i) suc/i 
i/iai Oj > (7j + i). 27iws i?„ /ias the same distribution as 1 plus the number of descents 
in a random permutation of [n] . 

Proof. Consider E n z Rn . Applying Lemma 14.11 k times and collecting the 
terms involving the expectation of the same expression we see that it is of the form 

1 k 

(4.2) E n z R " = — — V c k>m (z)E n _ k (z R ^(m + l) u ^), 

n ■ ■ ■ (n — k + 1) 

(for < k < n) with certain coefficient functions c k>m = c kim (z). Apply Lemma |4. II 
again to each of the expectations in the sum to get 

E n _ fc (z^-*(m + l)^-*) = (m + 1)0 f ~ 1} E M (z R ^-> (m + l) u ^) 

n — k 

+ !H±lE M (z R ^( m + 2)^—i). 

n — k 

Putting that back in, we see that (n ■ ■ ■ (n — k))gR n {z) is equal to 

k 

c fc , m (m + l)(z- l)E„_^ 1 ( 2 fl "- fc - 1 (™ + l) C/ "-' t - 1 ) 



m=0 



c fc , m (m + l)E n _ fe _! {z Rn -"-^{m + 2) u - ' ' i 

m=0 
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By rearranging the terms this is 



Cfc,m(m +!)(«- i)E n - k -. 1 (z B "-'- 1 (m+i) u "-*- 1 ) 



m— 1 

fc-1 



m— 

+ c M (fe + l)E„_ fc _ 1 (z ii - fc - 1 (fc + 2) c/ - fc - 1 ) 

k 

+ {c k , m (m + 1)0 - 1) + TOC fe , m _i} E n _ fe _! ( 2 r fl »-«- 1 (m + l) u »-*-i) 

m=l 

+ c fe . fe (fc + lJE^^fc*-'-'^ + 2) u »-*-i) 

fc+i 

= Y cj t4 .i, m E n _ fc _i(2r fl »-»- 1 (m + l) tr »-*- 1 ), 

where the coefficients Cfc !m satisfy the following recurrence: Cfc iJn = unless < 
to < fc, co,o — Ij an d for < m < k, 

Ck,m = mCk-l,m-l + ( m + 1)0 - l)Cfe— l,m- 

It follows by induction that c k<m {z) is a constant times (z — 1) m , so 

Ck,m — ^fe,mO 1) ; 

where ao,o = 1 and 

Jmak-i, m -i + (m + l)a fc _ M „, if < m < k, k > 1; 
1 0, otherwise. 

If we now let A; = n — 1 in (|4.2p and use Ei {z Rl (m + l)* 71 ) = z(to + 1) then we 



(4.3) o fc , 



n— 1 n— 1 



get 

(4.4) g Rn {z) = — V (to + l)cn_i, m = ^- V(m + l)a n _i. m (z - l)"" 1 " 
n! * — ' 7i! z — ' 

m— m— C 

The recurrence (|4.3|) is solved by 

. ffc + 1 

ak, m = ml < ,, 
{m + 1 J 

where { } is the number of partitions of the k element set into m non-empty 

subsets (recall the basic recurrence { m +i} = (m + l){ m+1 } + { m })- Hence, (|4.4|l 
becomes 

n-l , ^ 

fffl-(*) = iE( m+1 )M 1 1 (--i)- (m+1) 

n ' to + 1 

m=0 k J 



Z 



10 



PAWEL HITCZENKO AND SVANTE JANSON 



where in the last step we shifted the index by one and used {q} = for n > 1. 
This completes the proof since for a complex z (see (Kj Section 5.1.3]) 

^ \r - 1/ nJ. _ n\ ^ ml \m\ ^ ~ ^ ' D 

r=l N ' m=0 - 1 

Once the coefficients have been identified as Eulerian numbers we can use their 
known properties (see e.g. [DBi Chapter 10, pp. 150-154]) to obtain 

Corollary 4.3. The number of rows R n in a random permutation tableau of 
length n satisfies 

(4.5) ^-(n + l)/2^ 

V(rc + 1)/12 

In fact, a local limit theorem holds as well (see e.g. (CKSSj or [E]). 

Moreover, for the number of columns in a random permutation tableau, by 
definition, we have C n — n— R„ . Further, let D n denote the number of descents in 
a random permutation; thus R n = D n + 1 by Theorem l4.2l Hence C n — n — R n 
a I D n = A n , the number of ascents in a random permutation. By symmetry, 
A n = D n = R n — 1, and thus we obtain the following symmetry property: 



Corollary 4.4. The number C n of columns in a random permutation tableau 

satisfies C n = R n — 1. 
In particular, 

V(n + 1)/12 

The fact that R n — 1 = C„ = D n = A n follows also by the bijections described 
in [CoN] , 

Remark 4.5. The coefficients in (|4.4|) can be more explicitly written as 
(4.6) o B -i J n-r- l = (n-r-l)! Ji(j2-l)---(jr-(r-l)). 

l<Ji<-<3r<n-l 

This can be seen by putting them in a Pascal type triangle 

s \ 

a l,0 a l,l 

/ \ / \ 



dn— 1,0 On-1,1 ■ ■ ■ On-l,n-2 Cln— l,n-l 

and observing that by (|4.3p . a move down (either SW or SE) from a coefficient 
o-k.rn has weight m + 1. The value of a given coefficient at the bottom is obtained 
by summing over all possible paths leading to it from the root ao,o the products 
of weights corresponding to the moves along the path. Any path contributing to 
ffln-i,n-i-rj < r < n — 1 has exactly r SW moves and if they are from levels 
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ji, . . . ,j r (with level at the top and n — 1 at the bottom), then the path has 
weight 

(n - 1 - r)\(ji + l)j 2 (j 3 - 1) • • • OV - (r - 2)), 
and (|4.6[) follows by renumbering the terms. 

Expression (|4.6[) gives a direct way to get the moments; from (|4.4[) the proba- 
bility generating function of R n — 1 is 

1 n_1 

9B. n -i(z) = — V(n - r)a, i _ lj „_ r _ 1 (z - l) r , 

rj I * — » 



71! 

r=0 



z=i n 



r!(n — r) 

~Q"n — l,n— r— 1 



1 



and thus the rth factorial moment of R n — 1 is 

E„((i?„ - l)(R n - 2) • • • (i? n - r)) = ^gi^-xii 

r! (ri — r)\ v-^ , . 

= , E ii(Ja-l •••OV- (r-1 ■ 

l<ji<-<jV<n-l 

In particular, 

1 ™ 1 n — 1 

E„i?„ = 1 + E„(i?„ - 1) = 1 + - V j = 1 + 

n ^-^ 2 

and 

var(i?„) = var{R n - 1) = E„((i?„ - l)(i?« - 2)) + E„(i?„ - 1) - (E n (i?„ - l)) 2 

2 ■ / ■ \ n - 1 (n - l) 2 

(n — 1 m 2 4 

1<J1<J2<»1-1 

(n- 2)(3n-5) n-1 (n-l) 2 _n+l 
~ 12 + ~~ 2 4^ _ 12 ' 

which explains the normalization in (|4. 5|) . 



5. Probability generating function of the number of superfluous Is 

Let S n be the number of superfluous Is in a random tableau of length n. In 
|CoH| calculations based on (|2.1|) and (|2.2[) were used to show that 

(5.1) E n , n = ("- 1 K"- 2 ) , 

In this section we use the same approach to derive the expression for the probability 
generating function of S n , see Proposition 15.31 or equation (|5.5[) below. Although 
the form of the probability generating function looks rather unwieldy, it can be used 
to compute (in practice low order) moments of S n in a relatively straightforward, 
albeit a bit tedious way. We will illustrate it by deriving the exact expression for 
the variance of S n (see Proposition l5.5j) below. However, there does not seem to be 
an easy way to derive the central limit theorem for S n directly from the probability 
generating function. For this reason, in the next section we will rely on a different 
approach to establish the CLT for S n . 
We write in this section 



fc=i 
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where Vk is the increase in the number of superfluous Is when the length of a 
tableau is increased from k — 1 to k. As in the preceding section let Gk be the 
row number of the topmost 1 (provided that the fcth move is W and counting only 
unrestricted rows) and recall that it is set to be Uk-i + 1 if the A;th move is S. Then 
by the same argument as in Lemma 13.21 



Pfc(G fc =i|^fc_i) = i, j = l,...,t4- 1; 

and 

P k (G k = Uk-i + l\Tk-i) = ^- i . 

Moreover, the joint conditional distribution of (U k , Vk) given Tk-\ is 

(5.2) C((U k ,V k ) | T k -i) = (Gk+BmiUk.! - G k ),Bm(U k ^ - G k )), 

in the sense that if we further condition on Gk = m, then the distribution equals 
that of (m+X, X) with X ~ Bin([/fc_i — m) (thus the two occurences of Bin(f/fc-i — 
Gk) in (15. 2p signify the same random variable); further we interpret Bin(— 1) as 0. 

Lemma 5.1. We have, for all complex z and w, and k > 1: 
E k (z v -w u * | Tk-i) 

1 fzw + l^"- 1 , ( 1 \ fw\U k -i S 



l + w(z-l)\ 2 / V l + w(z-l) 



and 

1 



E kiz ^) = _|__ EjI _ 1 ( z ^( ZW+ ir-> ) 

^IT^T))^^ 1 ^ 1 )}- 



Proof. Using ()5.2|) we have 

E k (z v »w u > | ^-x) = E k (z mn ^- G ^w Gk+mn ^- G ^ | 

17* -l 

= J2 Ek{(zw) m ^ u ^- m ^w m I(G k = m) | 



+ it; x+i;fe - 1 P fe (Gfc = C/ fe _i + 1 | ^-i) 
E(|) ro E fe ((z^--)|^ l)+ .(|) 



m— 1 



Efw\ m / zw + L\ /w\ i 
(2) (—J 

m=l v 7 
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l + w(z-T)\ 2 J V l + w(z-l) 



1 - 



n (!) 



which is the first formula. 

To prove the second write 



= E k E k (z s "w u " | Tt-i) = E k (z s "^E k (z v "w u " \ F k -i)) 



tfc(2 "W 

and use the first part and the usual reduction by (|2.ip from P^ to Pfc_i. □ 

For £ > set bg — bg(z) — X^=o ^ so * na ^ we nave 
(5.3) zh + l = b l+1 and 1 + (z - l)b e = z e+1 . 

If we substitute bg for w in Lemma 15.11 and use (|5.3[) we obtain the following. 

Lemma 5.2. For < k < n - 1 and £ > 0, 

+ (l-^- 1 )E n _ fc _ 1 (z s -*-^y-*- 1 )). 

Proposition 5.3. FFe ftave 

_^ n — 1 

9s n (z) = —/S^c n ^i^{z)bg(z), 

£=0 

where the coefficients c mj g — c m ^{z) satisfy Co y o = 1 and /or to > i/ie recurrence 

c m +i,i = c m /(l - z~ ~ x )bn (z) + c m ^_iz _£ 6i?_i(z). 

In particular, c m ,g = unless < £ < m. 

PROOF. Skipping the dependence of bts on z for notational convenience, we 
show first that for < k < n — 1 we have 

k 

S n -k h U n - k 



(5.4) E„z s " = E„z s <" - „... (n 1 _ fc + 1) E C ^™ E « 



.fc« fl -*6i; 



m=0 



Indeed, assume inductively (|5.4|) for some k > (fc = being trivial). Apply the 
previous lemma to each of the terms E n _ k z^ n ~ k b 7 n to get 

1 k 

^ ' m=0 

+ (1 - z-^-^E^-fc-xz^-*-^^-*-^} . 
Separating the sums, rearranging the terms, and collecting the coefficients in front 
of E n -k-iz Sn ~ k - 1 blir^ k l we obtain that E n z Sn is equal to 

f _ J E (cfc, m -i6 m -i2- m + c^„6 m (l-z-™- 1 )) E„_ fc _iz s »-* - 1 6^- fc - 1 }. 

' m=0 

This completes the inductive proof of (|5.4|) with the coefficients Cfc iTn given by the 
specified recurrence. Choosing k = n — 1 in (I5.4[) and using the observation that 

Fi 7 Sl h Ul — h 

(because Si = and C/x = 1) completes the proof. □ 
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Remark 5.4. Just as in the proof of Theorem l4.21 the coefficients Ck, m can be 
put in a Pascal-type triangle 

co.o 

/ \ 



Cl.O 



Cl,l 



/ \ / 

Cn— 1,0 Cn-1,1 



\ 



\ 



Cn— l,7i— 2 Cn — l,n— 1 

This time a SE move from a coefficient c m .£ has weight z bi{z) and a SW 
move has weight (1 — z~ e ~ 1 )bi(z). As was in the case of the number of rows, the 
value of a given coefficient at the bottom is obtained by summing over all possible 
paths leading to it from the root co.o the products of weights corresponding to the 
moves along the path. This can be used to obtain an explicit expressions for the 
coefficients (c„_i iin ). Any path from co ; o to c n -i,n-i-r, < r < n — 1, has exactly 
r SW moves and n—l — r SE moves. The total weight of SE moves (they are from 

C,0, C.,1, . . . , C. )n _2-r) is 

n— 2 — r n— r— 2 

n z-j-\(z)=z-( n z r ) n b^z). 

j=0 j=0 

The SW moves may be from c.^ ± , . . . , c./ r , for some 0<£i<---<£ r <n — 1— r. 
A SW move from c.^. has weight (1 — z~ t j~ 1 )bt j (z). Therefore, 

n— i — 2 r 



C„-l,„-l-r(z) = Z ( 2 ) Yl ^ fe ( Z ) 

fe=0 



E 

0<e 1 <---<l r <n-l-r j=l 



This gives an explicit expression for the probability generating function of 5„, 
namely 



(5.5) 

where 
(5.6) 



9s n (z) 



1 



r=0 



r (-r) 



6 fc (z) U n , r (z), 



0<£i<---<^,.<n-l-r j=l 



4 )^(^)- 



(Note that all negative powers cancel in (|5.5[) since gs n ( z ) 1S a polynomial in 
z.) Although this expression looks quite complicated and it is not clear to us at the 
moment how to deduce the asymptotic normality of S n from it, some information 
can be extracted from it. We will illustrate it by deriving an exact expression for 
the variance of S n . 



(5.7) 



Proposition 5.5. For n > 2, the variance of S n satisfies 

(n-2)(2n 2 + lln-f) 



var(5„) = 



360 
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PROOF. We compute the second factorial moment of S 

1 d 2 / n ~ X 

(5.8) E n S n (S n - 1) = ( y~l c "-i,fcftfc 



2=1 



Notice that every path contributing to c n _i )T -, for < r < n — 3 has at least three 
SW moves and so its weight will have at least three factors of the form 1 — z - - 7 . 
Hence if it is differentiated twice and evaluated at z = 1 it will vanish. It follows 
that the sum on the right hand side of (|5.8[) reduces to the last three terms. For 
k = n — 3 we have 

c n -i,n-3b n -3 = j — I &fe J b n -zA n ^{z) = G(z)A n ^(z), 



where we have set 



\k=0 



-. n— 3 

G(z) := z-(V)TTfe fe ( z) 



So, the second derivative of G(z)A n ^{z) is 

G"(z)A n>2 (z) + 2G\z)A' na {z) + G{z)A'^{z). 

Since 

l<l<m<n-2 

A na {\) = A' n2 {l) = so we only need G(1)A^ 2 (1). Now, 

1 "~ 3 1 

(5.9) G(l) = _JJ(fe + l) = _ -. 

n! - L - L nm-1 
fc=0 v ; 

To compute A^ 2 (l), writing 

(1 - z-')(l - z- m ) = hx[z), h^(z)b m ^(z) = h 2 (z), 

we see that 

ti{{l)h 2 (l) + 2/ii(l)/i' 2 (l) + /ii(l)/i 2 '(l) - ^'(1)^2(1) - 2£ 2 m 2 . 

Therefore, 



n — 2 m 1 ri — 2 



m— 1 

= -^ n ( n -!)("- 2 X 5n - 11) (2n - l)(2n - 3), 
180 

and combining this with (|5 . 9|) we obtain 

Id 2 (n-2)(5n-ll)(2n-l)(2n-3) 

— T"7T (Cn-l,ra-30n-3) 

n\ dz z 



180 



We next handle 



-(3) ( n ~ 2 



— ;C„-i , n _i&„_i = j— ( TT 6fc ) 6„_i 

n! n\ \ - L - L / 
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The last expression is the probability generating function of a sum of independent 
random variables Wo, . . . , W n , where Wq = —n(n — l)/2 and for 1 < j < n, Wj is 
a discrete uniform random variable on {0, . . . , j — 1}. So, with Y n = 2j=o Wj we 
have 



If- h 

• , n Cri — l,n— l"n — 1 

n! dz z 



^ = EY n (Y n - 1) =J2^r(Wj) + {EY n f - EY n 



Now, 



so that 



i=i 



l<j<n, 



1 n 

E^ = -]T(j-l) 



4=1 



Furthermore, for 1 < j ' < n 



1 



j-i 



Therefore, 



fc=0 



n(n — 1) n(n — 1) 



j-iy _ (j-i)(j + i) 

12 



ra(ra- l)(2ra + 5) 



72 



Finally, putting the above together we get 



n! g?z 2 



Cn— l.n— l^n — ] 



n(n— l)(2n + 5) n 2 (n — l) 2 n(n— 1) 



2=1 



72 
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n(n- l)(9n 2 - 5n + 46) 
144 ' 



It remains to handle ^jc n _i. Jl _2^n-2- We write 

Cn-l,„-2&n-2 = /0)At,lO), 

where 

n-2 

/(z) = z-ev) j] Mz) 

and, according to (|5.6[) . 



n-2 



A„,i(z) = ^(l-z-^ 1 )fe £ (. 



Then the second derivative of c n _i.„_2&n-2 is 

/"(z)A„, 1 (z) + 2f(z)A' nA (z) + f(z)A>; iA (z). 
At z = 1 the first product vanishes. For the remaining two first notice that 

1 ^.f(,)=^ 1 )ff^i, 

J 



(n-1)! 



i=i 
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is a generating function of a legitimate distribution function. Therefore, /(l) = 
(n — 1)! and /'(l)/(n — 1)! is the expected value of a random variable represented 
by f(z)/(n — 1)!. Since this expected value is 

' j=l J k=0 

we obtain 

= (n-l)(n-2)(rc-l)! _ 
It remains to compute the first two derivatives of A n ,i(z) at z = 1. 

= E (^-^ - z_£ ) + fc-i'*"'" 1 ) . 

so that 

(5.10) < 1 (l) = E^ = n(n " 1) 6 (2n " 1) - 

£=1 

Also 



^=i 

Since b'^l) = (£ - 1)1/2 we get 

<i(D - E ((/ - iK 2 - ^ + D) - - "("-y- 1 ) . 

t=i 6 

Hence, 

^c„_ 1>7l _ 2 6„_ 2 | 2=i = 1(2/'(1K ;1 (1) + /(1K ;1 (1)) 

(n-l) 2 (n-2)(2n-l) (n-l)(2n-l) 
~~ 12 3 

(n- l)(2n- l)(n 2 - 3n + 6) 
~~ 12 ' 

Combining all of these calculations gives 

(n - 2)(5n - ll)(2n - l)(2n - 3) 



180 

(n- l)(2n- l)(n 2 - 3n + 6) 



12 

nfn- l)f9n 2 



n(n-l)(9n 2 -5n + 46) 



144 

(n-2)(n-3)(5n 2 -n-16) 



720 

By the same argument 

1 d 



- <? - — — ( h 

-n^n — l,n— l^n 

1 



-1 T t-„-l,7l-2"ra-2 



2 = 1 



EY n + ±-(f'(l)A nA (l) + f(l)A' nA (l)) 
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Using A n ,i(l) = 0, /(l) = (n - 1)!, (j^TTO]) . and the value of EF„ we get 
n{n-l) (n - l)(2n - 1) _ (n - l)(n - 2) 

which conforms to ()5.1|) . Finally, 

var(S„) = ES„(S„ - 1) - (ES n ) 2 + ES n 

_ (n-2)(n-3)(5n 2 -n-16) / (n - l)(n - 2) \ 2 (n-l)(n-2) 
_ 720 V 12 J + 12 

_ (n-2)(2n 2 + lln-l) 
~~ 360 ' 

which proves Proposition [53] □ 



6. Asymptotic normality of S n 

In this section we provide a self-contained proof of the following 
Theorem 6.1. isn^oo we have 
(6.1) 5» - n 2 /12 ^ 

Proof. As we mentioned, the form of the probability generating function of S n 
obtained in the previous section does not seem to be convenient to yield the central 
limit theorem. For this reason, we will rely on a bijective result of Steingrimsson 
and Williams [SWj . According to their result the number of superfluous Is in a 
permutation tableau of length n is equidistributed with the number of occurrences 
of a generalized patterns 31—2 in a random permutation of [n). (An occurrence 
of a generalized pattern 31—2 in a permutation a is a pair 1 < i < j such that 
(Tj_x > (Tj > <Ji.) To analyze that quantity, it will be convenient to think of a 
random permutation as generated from a sample X\, . . . , X n i.i.d. random variables 
each being uniform on [0, 1] (the permutation is obtained by reading off the ranks 
of Xi, . . . , X n ). If for 2 < i < j < n we let Iij := Ix i - 1 >x j >x i then 

S n = 

2<i<j<n 

Notice that from this representation we immediately recover (|5.ip since 



ES n = E £ I l , = ( n 2 1 )p(X 1 >X 3 >X 2 )= 1 -f n 2 1 

2<i<j<n \ ' ^ 

Similarly, we can easily obtain the asymptotic value of the variance: we write 

var(5 n ) = cov ( i »i,ji>- r t2,j2)) 

ii<ji 

and note that if {i\ — 1, *i, Ji} H {«2 ~ 1,*2,J2} = then h 1 .j 1 and ii 2j2 are in- 
dependent and so their covariance vanishes. In the complementary case, the main 
contribution comes from the cases that contribute 8(n 3 ) terms to the sum. We 
obtain 
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var(S n ) ~ —\cov(I 2 , 3,12,4) + cov (^2,5, h,s) + cov(7 2 ,4, ^3,5) 

+ COv(/ 2 ,5,^3,4) + COv(/ 2 ,3,^4,5) + COv(7 2 ,4, h^fj , 

as all other cases contribute 0(n 2 ) terms to the sum. We calculate: 

EJ a ,3 n h,4 = P(*i > X 3 > X 2 , X X >X 4 > X 2 ) = 2P(Xi > X 3 > X 4 > X 2 ) = ^ 



so that 



Similarly, 



1 fl^ 2 



C0V(I 2 , 3 ,I 2 ,4) = Y2 I 6 



E/ 2l5 n J 4 , 5 = P(*l > *5 > *2, *3 > *5 > X A ) 

= AP(X 1 >X 3 >X 5 >X 2 > X 4 ) = ^, 

E/2,4 n / 3 ,5 = P(Xi > x 4 > x 2 .x 2 >x 5 > x 3 ) = 

E/ 2l5 n i 3 , 4 = p(Xj > x 5 > x 2 , x 2 > x 4 > x 3 ) = 
e/ 2 , 3 n 14,8 = P(Xi >x 3 > x 2 , x 3 >x 5 > x 4 ) 

= 3P(Xj > X 3 > X 2 > X 5 > X 4 ) = i_, 

E/ 2 ,4 n / 4 ,5 = P(^l > *4 > *2,*3 > ^5 > X 4 ) = -L. 



Hence, 

(6.2) var(S„) 



111 1116 
l^ + ^ + l^ + 120 + 40 + 40~ 36/ 180" 



Remark 6.2. The exact value of the variance could be obtained by computing 
the other terms. 

Finally, to establish (|6.ip we will rely on results presented in [J] and [JLR1 
Section 6.1]; see |Bo] for a closely related theorem proved by the same method. See 
also [E] for a related simple proof of the asymptotic normality of the number of 
descents, cf. Corollary 14.31 above. We let A = A n := {a = : 2 < i < j < n}. 
Then S n = J^aeA^a- R- ecau that a dependency graph L for {I a : a S A} is 
any graph whose vertex set is A and which has the property that if V\ , V 2 are two 
disjoint subsets of A such that L has no edges with one endpoint in V\ and the other 
in V 2 then the families {I a : a G Vi} and {I a : a € V 2 } are mutually independent. 
For our purposes it is enough to consider L defined by the following rule: we put an 
edge between ai = and a 2 = (i 2 , j 2 ) iff {h - l,n,ji} H {i 2 - l,i 2 ,j 2 } ^ 0. 

If a\, . . . , a r € A then the closed neighborhood of {ai, . . . , a r } in L is defined by 

r 

N L (a!,...,a r ) = {J{/3 € A : (3 = or oti/3 € E(L)}, 

i=l 

where E(L) denotes the edge set of L. Note that in our case for every fixed r > 2 
|iV"i(ai, . . . ,a r _i)| = O r (n), where O r ( ■ ) means that the constant may depend 
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on r. Hence, trivially 



Ef/Jia!,. ■ • ,I aT ^) = O r (n). 



aGAfi(ai,...,a r _i) 

Since, cf. (|5.ip for the exact value, 

J2 E/ Q < |A| - 0(n 2 ), 

by |JLR[ Lemma 6.17] we conclude that the rth cumulant of S n (defined by 

K>r{S n ) = i~ k -$tk l°g0S„(Q)) where (j>x{t) — gx(e lt ) is the characteristic function) 
satisfies 

\nr(S n )\ = O r (n 2 ■ n r - 1 ) = O r (n r+1 ). 

Hence, for r > 3 



S n - n 2 /12 



= O r (n r+1 nT* r ) =o(l), 



\ v / n 3 7l80 

as n — > co. Since Ki(X) = EX and KaPO = var(X) we have 

and the theorem follows by the cumulant convergence theorem (see e.g. [JLR1 
Corollary 6.15]). □ 

Remark 6.3. Our results can be used to draw conclusions about some other 
parameters. For example, let Y n be the number of Is in the random permutation 
tableaux T n . Although we have not computed an explicit formula for the distribu- 
tion of Y n , we can easily obtain its asymptotic distribution. Namely, as n — > oo we 
have 

This follows immediately from Theorem 16. II upon noting that 

Y n = S n + C n = S n + 0(n). 

It seems straightforward to prove a similar central limit theorem for Z n , the 
number of Os in the random permutation tableaux T„, using further bijective results 
by |SW| ; the main difference is that we need to consider several generalized patterns 
simultaneously and the joint distribution of their numbers of occurrences. We leave 
this to the reader. 

References 

[BaS] E. Babson, E. Steingn'msson, Generalized permutation patterns and a classification of 
the Mahonian statistics, Sem. Lothar. Combin., B44b (2000) 18 pp. 

[Bo] M. Bona, The copies of any permutation pattern are asymptotically normal, preprint 

2007. arXiv:math. CO/0712. 2792vl. 

[BrCPRE] R. Brak, S. Corteel, R. Parviainen, A. Rechnitzer, J. Essam, A combinatorial deriva- 
tion of the PASEP stationary state, Electron. J. Combin., 13 (2006) R108. 

[BrE] R. Brak, J. Essam, Asymmetric exclusion model and weighted lattice paths, J. Phys. 
A 37, no. 14 (2004) 4183-4217. 

[Bu] A. Burstein, On some properties of permutation tableaux, Ann. Combin. 11 (2008) 

355-368. 



ASYMPTOTIC NORMALITY OF STATISTICS ON PERMUTATION TABLEAUX 



21 



[CKSS] L. Carlitz, D. S. Kurtz, R. Scoville, O. P. Stackelbcrg, Asymptotic properties of Euler- 

ian numbers, Z. Wahrsch. Verw. Gebiete 23 (1972) 47-54. 
[CI] A. Clacsson, Generalized pattern avoidance, European J. Combin. 22 (2001) 961-971. 

[C1M] A. Claesson and T. Mansour, Enumerating permutations avoiding a pair of Babson- 

Steingrimsson patterns, Ars Combinatoria, 77 (2005) 17-31. 
[Co] S. Corteel, Crossings and alignments of permutations, Adv. in Appl. Math. 38 (2007) 

149-163. 

[CoH] S. Cortccl and P. Hitczcnko, Expected values of statistics on permutation tableaux, 
in 2007 Conference on Analysis of Algorithms, Discrete Mathematics and Theoretical 
Computer Science AH, 325-340, 2007. 

[CoN] S. Corteel and P. Nadeau, Permutation tableaux and permutation descents, Europ. J. 
Combin. 30 (2009) 295-310. 

[CoW] S. Corteel, L. Williams, Tableaux combinatorics for the asymmetric exclusion process, 
Adv. in Appl. Math. 39(3) (2007) 293-310. 

[CoWl] S. Corteel, L. Williams, A Markov chain on permutations which projects to the PASEP, 
Int. Math. Res. Not. IMRN (2007) article ID rnm055. 

[DB] F. N. David and D. E. Barton, Combinatorial Chance, London, Griffin, 1962. 

[DDM] B. Derrida, E. Domany, D. Mukamcl, An exact solution of a one- dimensional asym- 
metric exclusion model with open boundaries, J. Stat. Phys. 69 (1992) 667-687. 

[DEHP] B. Derrida, M. Evans, V. Hakim, V. Pasquier, Exact solution of a ID asymmetric 
exclusion model using a matrix formulation, J. Phys. A: Math. Gen. 26 (1993) 1493- 
1517. 

[DS] E. Duchi, G. Schacffer, A combinatorial approach to jumping particles, J. Combin. 

Theory Ser. A 110 (2005) 1-29. 

[E] C.-G. Esseen, On the application of the theory of probability to two combinatorial 
problems involving permutations. Proceedings of the seventh conference on probability 
theory (Brasov, 1982), 137-147, VNU Sci. Press, Utrecht, 1985. 

[F] W. Feller, An Introduction to Probability Theory and Its Applications. Vol. I. Third 
edition, John Wiley & Sons, Inc., New York-London-Sydney, 1968. 

[J] S. Janson, Normal convergence by higher semi-invariants with applications to sums of 

dependent random variables and random graphs. Ann. Probab. 16 (1988), 305—312. 

[JLR] S. Janson, T. Luczak, and A. Rucinski, Random Graphs. Wiley, New York, 2000. 

[K] D. N. Knuth, The Art of Computer Programming, vol. 3. Addison- Wesley, 1996. 

[P] A. Postnikov, Total positivity, Grassmannians, and networks, preprint (2006). 

arXiv:math. CO/0609764. 

[R] J. Riordan, An Introduction to Combinatorial Analysis. Wiley Publications in Mathe- 

matical Statistics John Wiley & Sons, Inc., New York; Chapman & Hall, Ltd., London, 
1958. 

[S] T. Sasamoto, One- dimensional partially asymmetric simple exclusion process with 

open boundaries: orthogonal polynomials approach, J. Phys. A: Math. Gen. 32 (1999) 
7109-7131. 

[SW] E. Steingrfmsson, L. Williams, Permutation tableaux and permutation patterns, J. 

Combin. Theory, Ser. A 114 (2007) 211-234. 
[USW] M. Uchiyama, T. Sasamoto, M. Wadati, Asymmetric simple exclusion process with 

open boundaries and Askey- Wilson polynomials, J. Phys. A: Math. Gen. 37 (2004) 

4985-5002. 

[W] L. Williams, Enumeration of totally positive Grassmann cells, Adv. Math, 190 (2005) 

319-342. 

Department of Mathematics, Drexel University, Philadelphia, Pennsylvania 19104 
E-mail address: phitczenko<Smath.drexel.edu 

Department of Mathematics, Uppsala University, Box 480, SE-751 06 Uppsala, Swe- 



E-mail address: svante.jansonamath.uu.se 



